Nucleosynthesis in Neutrino-Driven 
Super novae 



C. Fr6hlich a W. R. Hix 13 ' 1 G. Martinez-Pinedo c 
M. Liebendorfer a F.-K. Thielemann a E. Bravo d K. Langanke c 

N. T. Zinner 6 

a Departement fur Physik & Astronomic, Universitat Basel CH-4056 Basel, 

Switzerland 

h Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831-6354, 

USA 

c Gesellschaft fur Schwerionenforschung, D-64291 Darmstadt, Germany 

Departament de Fisica i Enginyeria Nuclear, Universitat Politecnica de 
Catalunya, E-08034 Barcelona, Spain 

e Institute for Physics and Astronomy, University of Arhus, DK-8000 Arhus C, 

Denmark 



Abstract 

Core collapse supernovae are the leading actor in the story of the cosmic origin of 
the chemical elements. Existing models, which generally assume spherical symmetry 
and parameterize the explosion, have been able to broadly replicate the observed 
elemental pattern. However, inclusion of neutrino interactions produces noticeable 
improvement in the composition of the ejecta when compared to observations. Neu- 
trino interactions may also provide a supernova source for light p-process nuclei. 
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1 Modeling Core Collapse Supernova Nucleosynthesis 

The complexity of neutrino transport and the frequent failure of self-consistent 
models for core collapse supernovae to produce explosions have generally 
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divorced modeling of core collapse supernova nucleosynthesis from model- 
ing of the central engine. Nucleosynthesis simulations commonly replace the 
central engine of the supernova with a parameterized kinetic energy piston 
f Wooslev and Weaver . 19951 : [Rauscher et al. . 20021: Limongi and Chief 1 2003) 
or a thermal energy bomb ( Thielemann et al. . Il996l : iNagataki et al. . 19981 ). 
The energy of this blast wave, together with the placement of the mass cut 
that demarcates ejecta from matter that is assumed to fall back onto the 
neutron star, are tuned to recover the desired explosion energy and ejected 
56 Ni mass. These two methods are largely comp atible with the largest dif fer- 



ences coming in the inner regions of the ejecta ((Aufderheide et all 119911 ) . It 



is this inner region, where much of the iron, nickel and neighboring nuclei are 
produced, that is also most affected by the details of the explosion mecha- 
nism, including the effects of interactions between nuclei and the tremendous 
neutrino flux. 



While the importance of neutrino interactions is manifest in the name of the v- 
process and well documented for the r-process, neutrinos potentially impact all 
stages of supernova nucleosynthesis. During explosive nucleosynthesis, in the 
inner layers of the ejecta, where iron group nuclei result from a-rich freeze- 
out, interactions with neutrinos alter the neutronization, changing the ulti- 
mate composition. Galactic chemical evolution calculations and the relative 
neutron-poverty of terrestrial iron and neighboring elements place strong limits 
on the amount of neutronized material that may be ej e cted into the interstel- 
lar medium by core collapse supernovae ( Trimble! . ll99ll ). lHofTman et al.l (jl99fil) 
placed a limit of 10~ 4 M Q on the typical amount of neutron-rich (Y e < 0.47) 
ejecta allowed from a core collapse supernova. Those multi-dim ensional simu- 
lations of the central engine t hat produce explosions, (see, e.g.. iHerant et all 
1994 : Janka and Miiller . 1996| ) predict the ejection of much larger quantities of 



neutron-rich iron group elements than this limit. In an effort to compensate, 
modelers have been forced to invoke the fallback of a considerable amount 
of matter onto the neutron star, occurring on a timescale longer than was 
simulated. One common p r operty exhibited by recent multi-group simu l ations 
f Liebendorfer et al. . 2001 : Rampp and Janka . 20021: Thompson et al. . 20031: 



Buras et al. . 20031 ) is a decrease in the neutronization of the inner layers of 



the ejecta due to these neutrino interactions. This is a feature that current pa- 
rameterized nucleosynthesis models can not replicate because they ignore the 
neutrino transport. While the decreased neutronization seen in multi-group 
transport models would reduce the need to invoke fallback, it also makes any 
fallback scenario more complicated, since the most neutron-rich material may 
no longer be the innermost. 
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2 The Nucleosynthesis Implications of Neutrino Interactions 



Because of the impact of the neutrinos on the nucleosynthesis, the nucle- 
osynthesis products from future explosion simulations (utilizing multi-group 
neutrino transport) will be qualitatively different from parameterized bomb 
or piston nu cleosynthesis mode l s. Thi s was demonstrated by exploratory cal- 
culations by iMcLaughlin et a~ The dominant processes are vjv and 
e 1 * 1 captures on shock dissociated free nucleons, though at later times (and 
in regions with cooler peak temperatures) the more poorly known vjv and 
captures on heavy nuclei may contribute significantly. In addition to their 
impact on the electron fraction, these interac tions, as well as neutral cu rrent 
inelastic neutrino scattering off these nuclei ((Bruenn and Haxtonl . 1199 1). are 
also important to the thermal balance, potentially affecting the a-richness 
of the ejecta, there by altering the abundance o f important nuclei like 44 Ti, 
57 Fe, 58 Ni and 60 Zn (jWooslev and Weaverl . fl995h . iTravaglio et al l (120031) used 
a parameterized neutrino luminosity and spectrum ( Janka and MullerT 19961) 
to drive a supernova explosion and a tracer particle approach with a large 
nuclear network to examine the nucleosynthesis. They found significant im- 
pact of nuclear electron capture since some ejected zones reached densities 
> 10 8 gcm -3 . However, since these authors ignored neutrino captures, their 
simulations tell only half of the story. 



We have examined the effects of both electron and neutrino capture s in the 



contex t of recent multi-group supernova simulations. These models (see lFrohlich et al. 



2005al for more deta ils) are based on ful l y gene r al relativistic , spher ically sym- 



metric simulations ( Liebendorfer et al. . 2001 ). Pruet et al. ( 2005| ) have per- 
formed similar simulations using tracer particles from two dimensional simu- 



lations fBu ras et all 120 03). In both cases, artificial adjustments to the sim- 



ulations were needed to remedy the failure of the underlying models of the 
central engine to produce explosions. Also in both cases, the simulations were 
mapped to simplified models as later times, because the neutrino tranport 
simulations could not be run to sufficiently late times. While both of these 
shortcomings need to be addressed, these simulations nonetheless reveal the 
significant impact that neutrino interactions have on the composition of the 
ejecta. 



We observe three distinct phases in the evolution of the electron fraction of the 
matter that will become the innermost ejecta as it collapses, passes through 
the stalled shock and is driven off by neutrino heating. During core collapse, 
the electron fraction in these lower density, silicon-rich, regions is little changed 
either by the electron capture that is deleptonizing denser regions or the rel- 
atively weak neutrino flux. However, the combination of the larger neutrino 
flux after core bounce and the burning of silicon to iron in the still infalling 
matter greatly enhances the neutrino capture rates. With most of the mat- 
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Ca Sc Ti V Cr Mn Fe Co Ni Cu Zn 



Fig. 1. Comp a rison o f elemental abundances for Ca to Zn between models by 
Frohlich et al, ( 2005ah (circles) and Thielemann et al. Jjquares) and obser- 

vation determinations for metal-poor dGratton and SnedenT 199 lh (upward pointing 
triangles) and extremely metal-poor jCavrel et al. ~ l2004h (downward pointing tri- 
angle) stars. 



ter still tied up in relatively inert heavy nuclei, the greater abundance of free 
protons over free neutrons allows antineutrino captures to dominate, lowering 
Y e . The passage of the matter through the stalled shock raises the tempera- 
ture, dissociating the nuclei, but the concomitant increase in density prevents 
the lifting of the electron degeneracy. As a result of the high electron chem- 
ical potential, the balance of electron and positron captures strongly favors 
lower electron fraction. However, the combination of neutrino and antineu- 
trino captures favors higher Y e because of the slight dominance of neutrinos 
over antineutrinos as well as slightly higher abundance of neutrons compared 
to protons in the fully dissociated, mildly neutron-rich matter. As a result, 
the electron fraction undergoes only mild excursions in this phase. Eventually, 
continued neutrino heating (or perhaps some other mechanism) is sufficient 
to reenergize the shock, in the process lifting the electron degeneracy in this 
innermost eject a. As a result, the rate of electron captures drops while the 
rate of positron captures increases causing Y e to rise. While the dominance 
of neutrino captures over antineutrino captures drops as the matter becomes 
neutron-poor, their sum continues to favor higher Y e . Eventually, the electron 
chemical potential drops below half the mass difference between the neutron 
and proton, allowing p ositron and neutrino captures to dominate electron and 
anttaeutrino captures ' bsktotot B 

. With both neutrino emission and 
absorption processes favoring a higher electron fraction, Y e rises markedly in 
this phase, reaching values as high as 0.55. 
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interactions (iFrohlich et al 



Isotopic 



abu ndance s including the effects of neutrino 
2005bl ) relative to solar abun dances (filled circles) 



compared with earlier predictions ( Thielemann et al. . 1996h (open circles), which 
neglected neutrino interactions. The effect of neutrino interactions is clearly seen 
for nuclei above A > 64 where enhanced abundances are obtained. The inner panel 
presents the elemental abundances obtained from the isotopic abundances. 

The global effect of this proton-rich ejecta is the replacement of previously 
documented overabundances of neutron rich iron peak nuclei (nea r the N=50 
closed shell) ()Wooslev and Weaver , 19951 Thielemann et al. , 1996f ) with a mix 
of 56 Ni and a-particles. Production of 58>62 Ni is suppressed while 45 Sc and 
49 Ti are enhanced. Elemental abundances of scandium, cobalt, copper and 
zinc are significantly closer to those observed (s ee Fig. 1). T he results are 
however sensitive to the details of the simulations. IPruet et al. found a signifi- 
cant sensitivity in the nuclear production to the expansion rate of the matter, 
which was a parameter in their late time extrapolation. In addition to the 
global effects on the neutronization and entropy of the matter, our simula- 
tions, which include neutrino and antineutrino capture rates on heavy nuclei 
( Zinner and Langankel . 2004 ). find that these reactions have direct impact on 
the abundances of species like 53 - 54 Fe, 55 > 56 > 57 Co, 59 Ni and 59 Cu at the 10-20% 
level. 



A mild rp-process is also seen in this proton- rich ejecta (see Fig. 2). |Pruet et al. 



who neglected neutrino captures in their network, find this effect limited to 
zinc (A=64) by /3-decay lifetimes of waiting-point nuclei that are longer than 
the expansion timescale. We have found that transformation of protons into 
neutrons by neutrino captures allows (n,p) reactions to take the place of (3- 
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decays allowing significant flow to much higher A. We term this the z/p pro- 
cess because of the essential role the neutrinos play in producing these light 
p-process nuclei. The quantity of p-process nuclei produced and upper mass 
limit of this production are quite sensitive to the strength of the neutrino 
interactions and therefore the details of the neutrino source as well as the 
proximity and duration of the neutrino exposure. 



3 Conclusion 



Our results, and those of iPruet et al. (2005). clearly illustrate the need to in- 



clude the full effect of the supernova neutrino flux on the nucleosynthesis if 
we are to accurately calculate the iron-peak nucleosynthesis from core collapse 
supernovae. The sensitivities displayed in these models point strongly to the 
need to couple simulations of core collapse supernova nucleosynthesis to mod- 
els of the explosion mechanism. Not only will this foster better understand 
the contribution of supernovae to our cosmic origins, but comparison of nu- 
cleosynthesis estimates with observations will also improve our understanding 
of the central engine. 
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